Simulated synchrotron and Inverse Compton emission from 

Pulsar Wind Nebulae 



Delia Volpi*, Luca Del Zanna*, Elena Amato 1 ^ and Niccolo Bucciantini** 

* Dipartimento di Astronomia e Scienza dello Spazio-Largo Fermi 2-50125 Firenze-Italy 
' INAF-Osservatorio Astrofisico di Arcetri-Largo Fermi 5-50125 Firenze-Italy 
"Astronomy Department-University of California at Berkeley-601 Campbell Hall-Berkeley-CA 

94720-3411 -USA 

Abstract. We present a complete set of diagnostic tools aimed at reproducing synthetic non-thermal (synchrotron and/or 
Inverse Compton, IC) emissivity, integrated flux energy, polarization and spectral index simulated maps in comparison to 
observations. The time dependent relativistic magnetohydrodynamic (RMHD) equations are solved with a shock capturing 
code together with the evolution of the maximum particles energy. Applications to Pulsar Wind Nebulae (PWNe) are shown. 
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Pulsar Wind Nebulae (e.g. Vela and Crab Nebula) are 
a class of Supernova Remnants (SNR), which originates 
from the interaction of the pulsar wind with the surround- 
ing ejecta. They are characterized by non-thermal radia- 
tion at all wavelengths, mostly synchrotron (from radio 
to X-ray bands) and Inverse Compton (gamma-ray band), 
due to the presence of high energy pairs embedded in a 
strong magnetic field. In the optical and X-ray ranges ob- 
servations from space (HST and Chandra) show typical 
axisymmetric features known as "jet-torus structure". 

The theoretical interpretation of the jet-torus morphol- 
ogy Hfl, ifloll is based on the idea of a stronger equa- 
torial energy flux. This creates the torus and an oblate 
termination shock (TS) with cusps at the poles. Post- 
shock flows first converge toward the equator and then 
are diverted along the symmetry axis by magnetic hoop- 
stresses to form the jets. Thanks to the development of 
RMHD codes (i.e. |5[]), we can now investigate this pic- 
ture in detail. 

Our scheme |6] evolves the ideal axisymmetric 
RMHD equations together with the local maximum 
particle energy (in unit of mc 2 ) &», taking in account of 
adiabatic and synchrotron losses: 
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where the first term represents the adiabatic losses, the 
second one represents the synchrotron losses averaged 

with pitch angles, (j^) ~ 4 ° 
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fluid proper density, m and e are respectively the particle 
mass and charge, c is the velocity of light, t is the time 
and B is the amplitude of the magnetic field in the 
comoving frame. 



We employ the ECHO (Eulerian Conservative High 
Order Scheme) code developed by Del Zanna's code |5] 
(recently extended to GRMHD 18]) assuming a purely 
toroidal magnetic field and poloidal velocity. A descrip- 
tion of the initial and boundary conditions can be found 
in |@] and 0. 

Emitting particles are continuously injected at TS with 
a distribution function which is a power-law in energy 
and isotropic in momentum: 
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where A w is proportional to the thermal pressure of the 
post-shock fluid (the electron thermal energy), Ejs is the 
particle energy at TS and a = 0.7 is the spectral index. 
The particles are then advected by the nebular mildly 
relativistic flow and the post-shock distribution function 
is obtained from conservation of particles 'number along 
streamlines as in |7|]. 

Under the assumption of quasi-stationarity and of neg- 
ligible synchrotron losses, the emission coefficient in ob- 
server's fixed frame is: 
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if Voo > V and elsewhere. D is the Doppler factor, n is 
the direction of the observer in the comoving frame, v is 
the optical or X-ray observation frequency, v«, is the cut- 
off frequency (used to obtain the synchrotron burn-off): 
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Integrating the emission coefficient along the line of 
sight one can obtain surface brightness, optical polariza- 
tion [4] and spectral index maps 0]. 



Simulated flux maps reproduce the two polar jets and 
the equatorial flows with the correct ranges of velocity. 
Brightness maps show, as expected, an equatorial torus, 
a system of rings, with brighter arches and a central knot 
due to Doppler boosting (FigJT}. The optical and X-ray 
spectral index maps (FigQ} agree with observations lfl2tl . 

We found that, by changing the width of the equato- 
rial striped region b and the magnetization of the wind 
a (with deffective > 0.1 to have supersonic polar jets), the 
emission pattern changes so that it is, in principle, possi- 
ble to use the observed morphology of the jet-torus struc- 
ture to infer the conditions in the wind and to reproduce 
different Pulsar Wind Nebulae. 

High energy gamma-ray emission (COMPTEL, 
EGRET, HEGRA and HESS observations) is primarily 
due to IC scattering of high energy electrons on tar- 
get photons from CMB, dust (far infra-red emission, 
FIR), starlight, and from the PWN itself (synchrotron 
emission, SSC). 

We consider here two different post-shock distribution 
functions QJJ] : one describing the primordial radio emit- 
ting electrons (born at the supernova outburst) and the 
other the high energy tail accelerated at the TS. Our aim 
is to reproduce the observed spectrum in Crab. 

The low energy distribution function is: 

f r (e) = Al .e-( 2 * +1 ) e - £ / £ * (5) 

with A r a constant chosen to fit data at lGhz, a, = 0.26 
the spectral index and the radio particle energy cut-off. 
The high energy distribution function at TS is: 

f T s(£ TS ) = A W (£ TS + eo )-(2«w + l) e - £T s/ £ c (6) 

where A w is a constant proportional to the thermal pres- 
sure of the post-shock fluid, £o is the lowest wind parti- 
cle energy chosen to fit radio-optical data, a w = 0.7 is 
the spectral index and e c is the gamma particle energy 
cut-off at MeV frequencies. 

The incident photon density per unit of frequency is 
obtained either from a black body formula (CMB or 
FIR target) or from integrated synchrotron flux (SSC) 
assumed to be homogeneous and isotropic. 

Integrating distribution function and power per unit of 
frequency over the particle energy and incident photon 
frequencies for every energy regime [2], energy spectra 
and brightness maps are produced. 

The emission recipes has been successfully tested by 
comparing the broadband spectrum obtained with the 
Kennel and Coroniti spherical flow |9B, flJJ] against the 
Crab Nebula observations (fig. |2). However, when ap- 
plied to the flow structure from RMHD simulations, re- 
sults at synchrotron high energies are not satisfactory yet 
and further tuning of the emission model would be re- 
quired to fit all data (Figj2j. Gamma-ray surface bright- 



ness maps are shown in Figf3]for various sources of target 
photons. The strongest IC emission is due to SSC. 

The present work confirms jet-launching mechanism 
due to magnetic hoop stresses with the best agreement 
between simulations and Crab Nebula data given by the 
wind magnetization parameter deffective ~ 0.02. 

Further work is required to refine our model in order to 
compare with future Glast observations (paper in prepa- 
ration). 

This complete set for calculating simulated syn- 
chrotron emission, polarization and spectral index maps 
accounting for synchrotron losses can be used for other 
classes of objects (eg. AGN jets), in any scheme for 
RMHD (e.g. non-conservative, in full 3-D). 
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FIGURE 1. On the left: runA (b = 10 and a = 0.025) simulated brightness map in X-ray band, in logarithmic scale and 
normalized to the maximum value. On the right: runA simulated spectral index map in X-ray band. Distances from the central 
pulsar are reported on the axes, expressed in light year (ly) units. 
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FIGURE 2. On the left: Kennel and Coroniti 2-D spectra. On the right: runA spectra. Crab observed data are also plotted. 






FIGURE 3. Simulated brightness maps (in units of erg/ (sr cm 2 s Hz) in logarithmic scale) at 350GeV. On the left: IC from CMB 
target. In the middle: IC from FIR target. On the right: IC from synchrotron target. 



